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Pairing rotations and pairing vibrations are collective modes associated with a field, the 
pair field, which changes the number of particles by two. Consequently, they can be studied 
at profit with the help of two-particle transfer reactions on superfluid and in normal nuclei, 
respectively. The advent of exotic beams has opened, for the first time, the possibility to 
carry out such studies in medium heavy nuclei, within the same isotopic chain. In the case 
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studied in the present paper that of the Sn-isotopes (essentially from closed {Z - N - 50) 
to closed (Z - 50,N = 82) shells). The static and dynamic off-diagonal, long range order 
phase coherence in gauge space displayed by pairing rotations and vibrations respectively, 
leads to coherent states which behave almost classically. Consequently, these modes are 
amenable to an accurate nuclear structure description in terms of simple models containing 
the right physics, in particular BCS plus QRPA and HF mean field plus RPA respectively. 
The associated two- nucleon transfer spectroscopic amplitudes predicted by such model cal- 
culations can thus be viewed as essentially "exact". This fact, together with the availability 
of optical potentials for the different real and virtual channels involved in the reactions con- 
sidered, namely '^"'"^Sn-i-p, ^'*'^Sn+d and ^Sn+t, allows for the calculation of the associated 
absolute cross sections without, arguably, free parameters. The numerical predictions of the 
absolute differential cross sections, obtained making use of the above mentioned nuclear 
structure and optical potential inputs, within the framework of second order DWBA, taking 
into account simultaneous, successive and non-orthogonality contributions provide, within 
experimental errors in general, and below 10% uncertainty in particular, an overall account 
of the experimental findings for all of the measured ^'*''^Sn{p, t)'^Sn{gs) reactions, for which 
absolute cross sections have been reported to date. 
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I. INTRODUCTION 



Customarily, the fingerprint of shell closure in nuclei is associated with a sharp, step function- 
like distinction between occupied and empty single-particle states, in correspondence with magic 
numbers (l\ (for a recent example concerning '^^Sn, see [2J and [3J). At variance with the case 
of infinite, fermionic systems, in which there essentially exists a continuum of states at the Fermi 
energy, in finite-many-body (FMB) systems, like e.g. the atomic nucleus, sp is not, in principle, 
well defined, at least not in closed shell nuclei. This is because a sizable energy gap is observed 
between the last occupied (j<) and the first unoccupied orbitals. This is also the case for other 
FMB fermionic systems, e.g. Ceo fuUerene, which displays a relatively large HOMO-LUMO gap, 
of the order of 1.6 eV as compared to Sf ~ 15 eV H. In such cases, one can use as a working 
definition of O, (ej^ - SjJ/2 ((shomo - ^lumo)!'^ in the case of Ceo). Away from closed shells, 
medium-heavy nuclei become, as a rule, superfluid, the distinction between occupied and empty 
states being blurred around ep. The Fermi energy is, in such situation, well defined and equal 
to the energy for which the occupancy probability attains the value of one half. In keeping with 
this result, in the case of closed shell nuclei, sf can be properly defined as the minimum of the 
dispersion relation associated with pair addition and pair removal modes. 



II. PAIR-SPIN FORMALISM 



The mixing taking place in superfluid nuclei between particle and holes is economically em- 
bodied in the Bogoliubov-Valatin [[6l quasiparticle transformation 

ay = UyGy - VyOl = U'^Q^ - V'^q'^, (Ie) 

al = Ulal- V*ay = U[,a'J - Vy^, (lb) 

where 

Uy = uy* , K = v>"'^ (2) 

are the BCS occupation amplitudes, U'^ and V[, being real quantities, while a'^ = Qi(f>)alQ~^((p) = 
e~"^al and a'^ = Q{(t))ayQ'^ {(p) = e"^av are creation and annihilation operators referred to the 
intrinsic system of reference in gauge space (i.e. body fixed BCS deformed state), = exp{-iN(p) 
inducing a rotation of the angle in this (2-D) space (gauge transformation), being the number 
of particle operator (see Apps. A, B and C). The states |v) and |v), connected by the time reversal 
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operator, have the same energy (Kramers degeneracy). In keeping with the fact that (1) is a unitary 
transformation, 

u^ui + KK = + = 1 (3) 

The transformation (1) provides the rotation in Hilbert space of the creation and annihilation 
fermion operators, which diagonalizes the mean field, BCS pairing Hamiltonian (angle 6y, see 
Figs. A. 1(b) and (c)), 

{Hp - AN)mf = ^y^y - <^«o(^'^ + + <^«o' = ^^^^ + 

v>0 v>0 

where, 6y = Sy- A and 

Eys = Y^ey - Ey) + — , (5) 

y>0 

is the ground state energy (see Appendix |A]), while = . Li the above expression we find the 
particle number operators 

Ny = alay + alay, (6) 



the pair creation and annihilation operators. 

Pi = alat, Py = GyGy, ^ ) 

the complex (0 indicating the gauge angle) condensed (superfluid) Cooper field (see Appendix |B]) 
ao = (BCS \P\BCS ) = {BCS \P^\BCS )* = e~^''^a'^ = e'^'^ ^vK = J] 

v>0 v>0 

the quasiparticle energy 

Ey = {e^^+Ay\ (9) 
the absolute value (modulus) of the pairing gap 

A' = Ga'o = GY UX = Ge^"^ J]^*^^' = ^^^^ 

v>0 y>0 

and the quasiparticle number operators, 

Ny = alay + a'lay. (11) 

For the case of '^"Sn, the quantities 6y, Ey, U'^, V[ and U'yVy introduced above are given in 
Tables 1 and 2. 



II Pair-spin formalism 



5 



The ground state of the system referred to the laboratory system of reference 'K as well as to 
the intrinsic (body-fixed) frame 7C' (see Appendices A and B), can be written as 



|5C5(0)k = n «vav|0) = f]([/v + KaI4)|0) 



v>0 



v>0 



v>0 



v>0 



= |5C5(0 = O)V, 



(12) 



where Q. is the pair degeneracy of the single-particle subspace, leading to an overall (trivial) phase. 

The central feature of the \BCS) wavefunction, that is condensation of largely overlapping 
Cooper pairs, is captured by the pairspin (quasispin) formulation of superconductivity (see (81 19]1 
and refs. therein), which ascribes values (see Fig. [T]^a) and Appendix [A]) 



to empty states, and 



sAv) = Sy{v) = 0, s,(v) = -1/2, 



*.v(v) = sJv) = 0, s,(v) = +1/2, 



(13) 



(14) 



to occupied states of the non-interacting (normal) system. The symmetry axis in pairspin space 
(the z-axis) is referred to as the gauge axis (see Appendix \A\ Fig. A.l[ a)). A superfluid system 
(see Fig. 1 ;b), see also Fig. 



A.3 



Appendix A) is characterized by a collective pairspin S j_ = 
{Sj(,Sy}, which points in a direction perpendicular to z, associated with the azimuthal angle 2 0. 
This direction defines an intrinsic reference frame 7C' (see Fig. Bl, Appendix B), in which the 
components of the average total pairspin (5 ) = 2v>o('?(^)) the mean field ground state, take the 
values (see Fig. A.l ^b) and Appendix [A] see also (B28a) and B29a) 

{Sy')=0, 

v>0 

<5,) = <s.) = ^^(y7-c/;2). 



(15a) 
(15b) 

(15c) 



v>0 



Thus, 5x gives a measure of the mixing of empty and occupied states of the BCS solution 
(see Eqs. ( 13|14 ). This is tantamount to saying that 5x defines a privileged orientation in gauge 



space (see Appendix |A]), in keeping with the fact that ( [T2| ) is a wavepacket in the number of 
pairs of particles. An emergent property of the associated symmetry breaking phase transition is 
generalized rigidity in gauge space (see Appendix pi). That is, the system can be set into rotation 
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(or change its rotational frequency) in gauge space in terms of two-particle transfer reactions. The 
pairspin polarization may rotate collectively around the gauge axis, and the azimuthal angle Icp of 
5x (see Appendix [A] and Fig. A.l[ see also ( |B28a ) and (B29a)) is therefore a dynamical variable 



associated with pairing rotational bands. 

A. Order Parameter 

The modulus of the order parameter ([8]), that is, of the quantity a'^ ~ 2v>() UyV[, = A' /G, is, 
for medium heavy nuclei (A ^ 120) of the order of 1.4 MeV/G ~ 7, in keeping with the fact that 
G « 25 MeV/A (major j-shell approximation, see [fTOl| Chapters 2 and 3 and refs. therein). In 
other words, roughly of the order of ten (v, v) Cooper pairs contribute to the nuclear condensate in 
superfluid nuclei. Consequently, large fluctuations are expected for a'^. 

These fluctuations are generated by the residual interaction acting between the quasiparticles 
(cf. flSl, cf. also [[Toll in particular Ch. 4 and Appendices I and J of this reference). In the harmonic 
(QRPA) approximation the two associated pair fields are: 1) {U^ - Vy)(Tl + Y^) leading essentially 
to a bound two-quasiparticle like state (pairing vibration mixed to yS- vibrations, in deformed nuclei, 
cf. [fTT|| and refs. therein) lying on top of the pairing gap, = alal being the two-quasiparticle 
creation quasiboson operator; 2) (t/J -I- V^){YI - F^) which sets the \BCS)-k' intrinsic state into 
rotation in gauge space, and whose fluctuations diverge in the long-wavelength limit, in just such 
a way that the resulting ground state 

lA^o) - r d0e'^»l5C5(0)k - (J] c{y)ala^^'''l^\Q), (16) 

^ v>0 

transforms irreducibly under gauge transformation 

&{<P) = e-'^^ (17) 



The states (16) are the members of a pairing rotational band build out of a condensation of 
N/l Cooper pairs each described by |0) = Zy>of(v)aJajJ|0) (see Appendix B in particular Eq. 



(B16),( BT7| ) and (B23)). It is of notice that the presence of rotational bands in the spectra of 



many-body systems is the fingerprint of deformation both in real (3D) and abstract (e.g. gauge) 
spaces (see Fig. [2} see also Table XI in ref. [fT2l '). 

Fluctuations of the pair field are, of course, already present in the normal (correlated) ground 
state |0) = \gs{AQ)) of closed shell nuclei (mass number Aq). That is in systems in which, while 
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Ofp = 0, the (dynamical) value of the order parameter, i.e. the zero-point fluctuation of qtq around 
its zero value, 



1^ 



adyn = {(a - a'of) = 1/2 (<0|P+P|0) + {0\PP^\0)) = 

< \ 
\{int{A^ - 1)\pW + \{int(Ao + 2)|/'^|0>P 

int J 

^ ^ ({gs(Ao - 2)\P\0) + {gs(A + 2)\P''\0)) , (18) 

displays finite values (see [fTTl [T3l - fT5l and refs. therein), directly related to the quantity EcorrIG, 
where Ecurr is the average correlation energy of the pair addition {\gs{AQ + 2) > state) and of the 
pair removal {\gs{AQ - 2) > state) modes. In other words, the binding energy of the two-nucleon 
Cooper pair moving on top of the Aq closed shell system, and of the two-hole Cooper pair moving 
in the Aq Fermi sea. 



In the case of ^^^Sn (Aq = 132), we obtain (see Tablelinland Fig. [3 



E,orr{A + 2) = 2|e^> - ^1 - W{A + 2) = 1 . IVMeV, (19a) 

EcorAA - 2) = 2|e^< - .11 - W{A - 2) = 2.14MeV, (19b) 

and thus Ecorr = 1-66 MeV. It is of notice that in the case in which the energies W ^ 0, the system 
becomes superfluid, the BCS A parameter coinciding with the minimum value of the dispersion 
relation shown in Fig. 3|:a). Consequently £,o,,/G « ^ [^g^^^ + ^%^^] ~ ^(8.9+13.6) « 11, 
a value which is not very difl^erent from that of A' /G associated with the superfluid nucleus '^"Sn. 



The quantity ( 18 1 can be also estimated from the ratio (A/G) between the (two-particle)-(two- 
hole) (pairing) vibration coupling strength (see e.g. IfTTI ) and the pairing coupling constant. In 
the two-level model [fTDl , this quantity is given by (A/G) = 2 Vo, Q being the pair degeneracy of 
the single-particle space in which the two nucleons (two nucleon holes) participating in the pair 
addition and removal modes are allowed to correlate. As a rule, these are the valence shells of the 
closed shell system Aq. 

In the case of the newly discovered '^^Sn, doubly magic nucleus, Q a: 30 and a^yn « (A/G) « 



2 V30 « 1 1. Making use of the actual values of A (see Table III) one obtains | -I- ^] ~ 9. 
Summing up, adynC^^Sn) > aC^^Sn). 

This result embodies the very difi'erence between gauge spontaneous symmetry breaking in 
atomic nuclei and in condensed matter. In a chunk of e.g. Pb, of which more than 50% is the 
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atom built on the isotope ^''^Pb, (or of Sn, of which none is the atom built on the isotope '^^Sn, 
the corresponding nucleus being highly unstable), at a temperature below the critical temperature 
Tc = 1 .2K{3.12K) at which the metal becomes superconducting but in the presence of a magnetic 
field stronger than the critical value He = 0.08 T (0.03 T), of the order of 10^ times the earth 
magnetic field, Cooper pairs break as soon as they are formed, leading to a hardly observable 
effect, in particular concerning the structure and stability of the crystal. 

On the other hand, in the case of their ground state and thus at absolute zero temperature - as it 
is the case for all natural occurring nuclear species on earth - the Cooper pairs associated with the 
normal (non superfluid) system displaying two nucleons above or two holes below closed shell, 
like e.g. IgsQ^^Fhug) > and |g5'(g2^Pbi24) > respectively, these fermions are strongly correlated, as 
evidenced by the large two-nucleon transfer cross sections with which they are excited (see e.g. 
[fT2ll and refs. therein). In keeping with the fact that a consistent fraction of this cross section 
arises, e.g. in the case of ^^^Fh{t, p)^^'^Fh(gs) by the transfer of two particles to levels below the 
Fermi energy of 2o**pb (ground state correlations). Cooper pair correlations blur dynamically the 
difference between occupied and empty single-particle states thought to be a trademark of closed 
shell systems. 

The same arguments presented above, can be used for jQ^Sn82, in which case the summed 



backwardsgoing amplitudes amount to 2/ Yi ~ 0.5, as reported in Table III Arguably, one can 
posit, that in the nuclear case it is not, or at least not only, the condensed (superfluid) state which 
is peculiar, but the normal state [|T6l. in which pair addition and pair removal Cooper are virtually 
poised in the ground state of the closed shell system Aq, ready to condense (i.e. W(Ao + 2), W(Ao - 
2) 0, see Fig. |3]). inducing a (fluctuating) alignment of pairspins perpendicular to the gauge 



(z-axis), and thus to a domain wall and associated generalized rigidity (see Fig. A.3[ ). This can 



also be seen from the pairspin states in normal systems (see Figs. 4 and 5), defined as, 

|0)v = (YaM\s,(v) = -1/2) + XrMs,(v) = 1/2)) 

Norm. 

^ (Wv)|2), + X„,„(v)|l)v), foT€y<eF, (20a) 



Norm. 



and 



|0)v = iYrem{v)\s,(v) = +1/2) + X„rfrf(v)| 5,(v) = -1/2)) 

Norm. 

iY,e,n\l)v + X,M{v)\2)y) , for ey > ep, (20b) 



Norm 
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where X, Y are the forwardsgoing and backwardsgoing RPA pair vibration amplitudes (see insets 
ofFigjB]) 

Xn(v;/3) I ^ (Vn;/2)A„C6) 

The quantity Q.y = {2jy + l)/2 is the pair degeneracy of orbital v, /3 = -2 and yS = +2 label 
the pair addition and pair removal modes respectively, while n {= 1,2,... ) numbers the solutions 
of the RPA dispersion relation in subsequent order of excitation energy (we deal here only with 
the n = 1, lowest energy pairing vibrations). Consequently, as soon as the single-particle field 
Hsp - field which acts on the nuclear pairspin along the gauge (z)-axis in a similar way in which 
a magnetic field acts in the case of a metallic superconductor (see App. A, discussion after Eq. 
(A27); see also discussion after Eq, (A35)) - is decreased, a fact that takes place moving away from 
magic numbers (in the case of e.g. ^^^Sn this implies reducing the single-particle gap \sj^ - Sj>\ 
from about ^ 5 MeV to few hundreds of keV, se Tables|Il]and III), nuclear Cooper pairs condense, 
the ground state of the corresponding nucleus becoming amenable to a BCS-type description, a 
fact which already takes place with the presence of two removal modes in Aq. 

It is of notice the presence of strong fluctuations in pairspin observed in Fig. 5, as compared 
to the smooth variations shown in Fig. 1, idealization of the situation representative of a high 



purity, metallic crystal (within this context see Fig. |A.3| of Appendix [A]). In fact, finite nuclei 
display orbitals which contribute very differently to pairing corrleations, in particular hot orbitals 
(see e.g. IfT^fTTlfT^ ), related to the diff"erent pair degeneracy Qy and with the relative amount of 
i'-component of the diff"erent jy(0) pure two-particle configurations. One can view such an imper- 
fect pairspin alignment as a limitation of its applicability to finite many-body systems. Conversely, 
one can interpret it as a reflection of the richness with which these systems in general, and finite 
nuclei in particular, embody symmetry breaking phase transitions. Namely, among other things, in 
terms of very non-conventional normal phases. Normal phases which display (virtually) traces of 
e.g. domain walls with varied degree of stability, dynamically violating the symmetry in question 
(gauge symmetry in the present case). These properties are precisely those which render the study 
of pairing correlations in nuclei central, in the quest for the mechanism which are at the basis of 
the stability of nuclear species, in particular along the drip lines [|T9l . 

In keeping with the fact that the P"^ and P are the basic operators entering both the pairing inter- 
action {Hp = -GP'^^P') and the pair mean field {-Ga'^iP^ + P) + G{a'^^), two-nucleon transfer can 
be viewed as the specific probe of pairing correlations in nuclei, in a similar way as Coulomb exci- 
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tation, inelastic scattering and y-decay are specific tools to probe (quadrupole) surface vibrations 
and rotations (see Fig. 2). In other words, specific information on Uq and adyn can be obtained 
through two-nucleon transfer reactions. 

Because the correlation length associated with the nuclear Cooper pairs is 



tivp 
2A'' 

tiVf 



K * 0)' 

(22) 



IE 

-i-i^ con- 
Cooper pair partners are correlated over distances considerably larger than nuclear dimensions 

(^ « 20 - 30 fm, as compared to i? ~ 5 fm). Consequently, from a nuclear structure point of view, 

Cooper pair transfer involves also regions in which the pairing interaction G{x), x representing 

e.g. the surface-surface distance between target and projectile (see e.g. [12011 Ch. Ill), vanishes, 

a situation already known in condensed matter in connection with the Josephson eff"ect [BTI - I^ . 

This result, together with the fact that the depth of the single-particle potential |Vol(* 50 MeV) is 

much larger than G(« 25/A MeV), implies that, exception made for value eff'ects, successive 

transfer induced by the mean field single-particle (Saxon- Woods-like) potential, is expected to 

be, as a rule, the largest contribution to the two-nucleon transfer cross sections [|24] - l35lll36l . The 

difficulties to absorb this simple result by nuclear structure practitioners partially stems from the 

fact that, neglecting reaction details, the two-particle transfer cross sections gs gs can be 

schematically written as 

l«ol' (t^o ^ 0), 
(T{gs gs)~ I (23) 

\adyn? (^0 = 0), 

corresponding to the square modules of matrix elements of and P, the associated spectroscopic 
amplitudes in the intrinsic coordinate system in gauge space 7C' being (see e.g. [fT2ll . and Appendix 

m 

B{vv- 7 = 0) ~ ^^{BCSiA + 2)\P'J\BCS{A))^,, (24a) 

where 

TC'{BCS(A + 2)\P'J\BCS(A))^, = U'MW'M + 2) (= -x'iBCS {A)\P',\BCS (A + 2)kO . (24b) 
Now, thinking in terms of the transfer of a nucleon at a time, e.g. in the case of the reaction 

(A + 2)+p^ Fi=A + \) + d^A + t (25) 



A Order Parameter 
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one can write, in connection with the first step {(A + 2) + p ^ F + d), 

ala,\BCS (A + 2))\p) - K(A + 2)al]BCS (F))\d), (26) 

and 

a',]a',\BCS(F))\d) - f/;(A)K(A + 2)\BCS (F))\d), (27) 
in connection with the second step (F + d ^ A + t), in keeping with the fact that 



U'al + VyUy, 



4 = f/>y - Vial, 

and that \BCS ) is the quasiparticle vacuum. The primed quantities are referred to the intrinsic, 
body-fixed frame (see App. A and B). 

Consequently, the associated (successive) two-nucleon transfer spectroscopic amplitude 

5(y2(0)) = UMWM + 2) (= cy) , (28) 

has the same dependence on the BCS occupation numbers as that displayed by the amplitudes 
associated with simultaneous transfer (order parameter) [fT2|, and with the c(v) amplitude entering 
in the Cooper pair wavefunction. This last result is closely related to the smooth behavior of the 
BCS occupation parameters with mass number (see Appendix |B]). Summing up, pair coherence is 
maintained both in successive as well as in simultaneous transfer. 

It is of notice that all the above results, which constitute the very essence of nuclear BCS, 
are not only inescapable, they are also almost tautological, at least for well bound nuclei. In 
fact, pairing in nuclei does not affect neither the single-particle energies Sy (see in any case Eq. 



(A27) and following discussion), nor the corresponding wavefunctions (pvir), but only the single- 
particle occupation probabilities. And this takes place in a small (A/sf « 5 x 10"^) region around 
the Fermi energy. In this region, and in keeping with the structure of the BCS wavefunction, 
which takes into account the variety of excitations of pairs of particles (v, v'), so as to produce 
the most efficient mixing of empty and occupied states leading to Cooper pairs, the only possible 
excitation mechanism of the nuclear superfluid, is that of breaking a Cooper pair, individual two- 
particle (vv) excitations being already taken into account in the BCS ground state. It is then 
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neither surprising that the amplitudes of the Cooper pair wavefunction are c(v) ~ U[,Vy, nor that 
the absolute cross section for Cooper pair transfer are proportional to (liv>o U'yV^)^. And least of 
all, it is not surprising the fact these amplitudes and two-nucleon transfer processes involve not 
only (nljm,nlj - m) configurations, but also {nljm,n'lj - m) as well as [(J{A) + niliji)j^ ®(J'(a) + 
nihil) jM- This is in keeping with the non-orthogonality of the single-particle wavefunctions 
describing the target and projectile in the different channels {a + A ^ f + F ^ b + B). Within this 
context, not only n and n' are possible, due to the fact that partners of a Cooper pair feel different 
mean fields ^nijm)^ ^^^^ because a general nuclear structure treatment of pairing, will 

include Cooper- like correlations associated with multipole pairing (see e.g. JTOl Sect. 5.3 and 
refs. therein), correlations which, in the present case, have not a dynamical origin (one works with 
Hp = -GP^P), but only a trivial kinematical one. 



Ill, REACTION MECHANISM 



In what follows we present the elements which enter the calculation of the absolute two-particle 
transfer differential cross section in terms of the reaction 

A + t^ B(=A + 2)+p, (29) 

in which A + 2 and A denotes the mass number of even nuclei in their ground state. In other words, 
one concentrates on L = transfer. The wavefunction of nucleus A -I- 2 is written as, 

¥a+2(^a, Tai, CTi, r^2, o-i) = ipAiU) ^1' ^^2, o-2)]o . (30) 

product of the wavefunction describing the ground state of the nucleus A, the corresponding rel- 
ative (intrinsic) 3A - 3 radial coordinates being denoted and of the wavefunction of two- 
correlated nucleons 

(l>&Ai,CTi,rAi, cri)]l = Clnm [(pt}^j,{rAi, (Ti)(pifj^irAi, (Ti)]l (31) 

nm 

the wavefunctions i^^];^^ (r, cr) describing the single-particle motion of a nucleon in a mean field 
potential, e.g. a Saxon-Woods potential. The spatial part of the two-neutron wavefunction in 
the triton can be written as, 0,(rpi,rp2) = p{rp\)p{rpi)p{rii), r^i and Vpi denoting the modulus 
of the relative coordinate of each of the two neutrons involved in the transfer process, measured 
with respect to the proton, while ru denotes the modulus of the relative coordinate of the two 
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neutrons in the triton. The functions p(r), as depicted in Fig. [6]^a), are generated with the p - n 
Tang-Hemdon interaction [137]| 

v{r) = -vq exp (-kir - r^)) r > (32) 
v{r) = 00 r < rc, (33) 

where k = 2.5fm"' and = 0.45fm denotes the radius of the hard core. The depth vq is adjusted so 
as to reproduce the binding energy of the triton and of the deuteron respectively. This hard-core 
potential is also used in the above expressions as the n - p interaction potential responsible for 
neutron transfer. 

The two-particle transfer differential cross section is written as 

da _ MiMf kf I _ ,2 

dQ {Ann^fki I ■ ^ ^ 

The amplitudes appearing in it describe the simultaneous, 



liji 0-10-2 



X v(rpi)(f>,(rpi, rp2);rM (I'm), (35a) 



successive 



T'sil = 2Yj Z Z fdr,Fdrp,drA2[(f>t^(rAu(rurA2,cr2)t/;^\rpB>^^^^ 



xMrpi)(pf/;l^^irA2) J dr'^,pdr'p^dr'A2G{rdF,r'aF) 

X <^.(r;j>f^-^%/r;,2)^Kr;2)<^.(r;i)0,(r;2);r;iVM), (35b) 



and non-orthogonality 



t(2) _ 9 
^ NO ~ 



Yj Yj I drdFdrpidrA2[((>f_^jf(rAucri,rA2,o-2)]o*X^pB(^pB)virpi) 

l.J, lf,jf,mfO-iO-2 J 

xMrpM;:l„^^(rA2) J dr'p.dr'^^dr'. 



dF 

X 0.(r;i)>J;i:„/r:,2)'^Xr;i)<^,(r;2);rlI^(rM), (35c) 

contributions to the transfer process. In these expressions, ip^^^^^^^^{rA\) are the wavefunctions de- 
scribing the intermediate states of the nucleus F = A + l, generated as solutions of a Saxon-Woods 
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potential, while (pd^pi) is the wavefunction describing the deuteron bound state (see Fig. |6|b)). 
We have chosen the so called post-post representation [|20]| . in which the p - n interaction appears 
twice in the successive amplitude. The Green function G{rdF,r'^f) propagates the intermediate 
channel d, F, and can be expanded in partial waves 

G{rdF, r^p) = z > V2Z + 1 [Y (r^^ )y (r^^)J^ . (36) 

J l^dFf dFf idF 

The functions fiik^F, r) and Piik^F, r) are the regular and the irregular solutions of a Schrodinger 
equation associated with a suitable optical potential and an energy equal to the kinetic energy in 
the intermediate state. In most cases of interest, the result is hardly altered if one uses the same 
energy of relative motion for all the intermediate states. This representative energy is calculated 
when both nuclei appearing in the intermediate state are in their ground states. The validity of 
this approximation can break down in particular cases. For example, in the case in which some 
relevant intermediate states are strongly off shell, in which case their contribution is significantly 
quenched. An interesting situation can develop when this situation becomes operative for all 
possible intermediate states, in which case they can only be virtually populated, thus emphasizing 
the role of simultaneous transfer. 

IV. THE ISOTOPIC CHAIN ^° 'Snso-^,^ Sngi 

A collective mode is characterized by: 1) an enhanced cross section of transition probability; 
2) a simple expression of its energy as a function of the quantum number characterizing the states 
connected by the transition. This quantum number is related to restoration of the symmetry vi- 
olation (static or dynamical, e.g. particle number in the case of pairing rotations and vibrations, 
angular momentum in the case of e.g. quadrupole rotations and vibrations). 

For example, in the case of a quadrupole rotational band of a 3D-deformed nucleus like, e.g. 
'^^Dy, 1) corresponds to the B{E2) transition probability, measured e.g. in the terms of single- 
particle Weisskopf units (of the order of 10^ in the example chosen), while 2) corresponds to 
Ej = {fp- I2J')I{I + 1), / being the angular momentum of the system (/ = 0, 2, 4, ...). In the case of 
pairing rotational bands 1) corresponds to the absolute value of the two-nucleon differential cross 
section, measured in terms of the average pure two-particle units ffT2l [TTl [TSlI (typical value of the 
enhancement factor being, in the case of Sn-isotopes, of the order of 10^), while 2) corresponds to 
= {fp-l2I){N - Nof, N being the number of particles associated with the condensate, namely 



A Pairing rotations 



15 



with the neutrons in the case of jgSn/v, while A^o is the mean number of neutrons representative 
of the particle number wavepacket describing the superfluid Sn-isotopes (A^o * 68, see discussion 
below). Of these two quantities, 1) is arguably the most representative one. This is because 
accidental degeneracies or residual interactions may modify the energies without much altering 
the long-range correlation of the coherent state. B8I This is also the reason why, in what follows, 
use is made of the single j-shell model to discuss the basic features of the pairing rotational modes. 

In this Section we present evidence of the accuracy with which the model of pairing rotations 
discussed in Section 2, together with the two-nucleon transfer reaction scheme summarized in the 
last section allows for an overall quantitative description of the absolute value of the two-nucleon 
transfer cross sections, when use is made of global optical parameters to describe the (three) elastic 



channels involved in the process. Consequently, the predictions given in Sect. IV B concerning 
the pairing vibrational spectrum expected in connection with the two unstable closed shell systems 
^^^Sn and the most exotic one ^""Sn can be considered potentially important and likely quantitative. 



A. Pairing rotations 

In Fig. |7] we display the value of the absolute differential cross sections associated with the 
reactions ^Sn(;?, t)^QSn{gs) for which absolute measurements have been reported in the litera- 
ture, in comparison with the experimental data [|^9] - |45l (see also [|46l HTII ). The corresponding 



integrated cross section are collected in Table IV In all cases the contribution of the successive 



process is the dominant one. Examples of two-nucleon spectroscopic amplitudes obtained from 
BCS calculations are displayed in Table 2 (U, V for ^^"Sn(/?, O'^^^Sn). They have been computed 
solving the gap and number equations with a monopole interaction acting on the bound orbitals, 
calculated as the eigenfunctions of a standard parametrized Saxon-Woods potential, and imposing 
that the gap reproduces the value obtained from the empirical odd-even mass differences for the 
various isotopes. The BCS spectroscopic amplitudes are in good agreement with those predicted 
by extended shell model calculation (see refs. [|43ll44ll and refs. therein). The optical parameters 
in the entrance, intermediate, and final channel where taken from refs. in ||39] - l44ll and from [|48ll 
for the deuteron channel. 

From the above results one can posit that theory provides an account of the experimental ab- 
solute differential cross section well within the experimental errors and, arguably, without free 
parameters. 
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Let us now concentrate our attention on the value and on the structure of the probability ampli- 

-> 

tude for two nucleon, at r and r' to belong to a Cooper pair, namely a'^ir, f') = 2v>o Cv(pv{r)(t)y(r') 



(see (16) and (28)). That is, the nuclear structure component of the two-particle transfer cross 
section amplitude. To clarify the physics at the basis of the BCS description of pairing rotational 
bands, we discuss two scenarios for the case of the '^*'Sn(p, 0^'^Sn reaction. In the first one, we 
consider all bound single-particle states, the cutoff energy being E cutoff = MeV. In the second 
case one sets Ecuwff = 60 MeV, discretizing the continuum inside a spherical box of 15 fm of 
radius. The BCS gap (A = 1.47 MeV; experimental value) and number (A^ = 70) equations lead 
to G = 0.18 MeV and /{ = - 6.72 MeV in the first case and G = 0.05 MeV and A = - 6.9 MeV 
in the second one. The associated Cooper pair probability distributions in r-space are essentially 
identical (see Fig. [8] and Tables |l]and|ll]). It is then not surprising that they lead to essentially 
the same absolute value of the two-particle transfer cross section associated with the reaction 

Let us now repeat the argument, but this time in terms of cr(gs gs) ~ (A/G)^, as it customary 
done since the first publication which introduced it B9ll . Because the pairing gap has been fixed 
to reproduce the experimental value (1.47 MeV), one obtains in the case of Ecutoff = OMeV 
(A/Gf ~ 70 and (A/Gf ~ 889 in the case of Ecutoff = 60MeV. This result emphasizes the 
problem of working with an expression which contains explicitly the pairing coupling constant. 

One could argue that such an objection could also be leveled off against the relation cr ~ loroP- 
Note however, that a (p, t) reaction would hardly feel the effect of contributions far removed from 
the Fermi surface A. This is in keeping with the fact that transfer to levels lying far away from A will 
be unfavorable due to Q-value effects. If one argues in terms of the relative distance r between 
target and projectile (r » Rq for continuum-like contributions; r < Rq for deeply bound-like 
contributions), the outcome is similar. In fact, for large distances the two-particle transfer form 
factor vanishes while at small distances the outgoing tritium will experience strong absorption (see 
Appendix |E]). 

In fact, considering only the contribution to a,, arising from the valence orbitals, that is, es- 
sentially those contributing to the "naked" vision of the Cooper pair wavef unctions, one obtains 
a'^ = 2.12 and a'^ = 2.08 respectively, and thus, a negligible squared relative difference between 
the two predicted cross sections, namely (0.04/2.1)^ ^ 2x 10"^. 

Summing up, because the pair condensed state can be viewed as a coherent state which behaves 
essentially classically when viewed in terms of its building block (Cooper pair) the description of 
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pairing rotational bands provided by the BCS model in terms of a coupling constant and an energy 
cutoff can be considered essentially "exact" when probed with two-nucleon transfer processes, 
reactions which filter the inaccuracies of each individual UyVv component, emphasizing the off- 
diagonal long range order provided by the phase coherence (cf. ref. [ISOlD . In fact, studying nuclear 
Cooper pair condensation in terms of e.g. single-nucleon transfer (see e.g. [l5Tll52l . the individual 
inaccuracies of the BCS occupation numbers cannot be averaged out. As a consequence, the 
overall agreement between theory and experiment is much poorer than that reflected by e.g. the 



results collected in Table IV The above arguments provide further evidence of why two-nucleon 
transfer is the specific probe of pairing superfluidity. 

In Fig. |9]a quantity closely related to the Sn-isotopes binding energy is reported as a function 
of the number of neutrons. Also displayed is the best parabolic fit to these energies, a quantity to 
be compared with 

Em = ^(N-No)\ (37) 

namely the energy associated with the members of the pairing rotational band. 

A simple estimate of the pairing rotational band moment of inertia is given by the single j-shell 
model (see e.g. 1101 App. H n^/2I = G/4 « 25/(4A^o) MeV). 

This estimate turns out to be rather accurate, even beyond expectation. Of notice that to the 



extent that one is discussing properties of a coherent state hke that described by (12), for which 



Hsp plays a secondary role (see discussion following Eq. (A27) in Appendix Kl) this is not a 



surprising results. As can be seen from Fig. [9[ the estimate (37 1 is rather accurate except close to 
=50 and A'^ = 82, in keeping with the fact that, as discussed before, the pairing deformed picture 
(oq 4^ 0) breaks down around closed shell (qtq = 0), where a vibrational regime (associated with 



the dynamic distortion a^y,,, see Eq. ( 18 1) is expected to be valid. 

Also reported in Fig. |9| are the integrated values of the measured absolute two-particle transfer 
cross sections. Naively, one would expect a marked constancy of these transitions, in keeping with 
the fact that the (pairing) rotational model implies a common intrinsic (deformed) state (Cooper 



pair condensate, see Eq. (12)). On the other hand, due to the fact that the number of Cooper 
pairs contributing to the pairing distortion «[, is rather small (less than 10), one expect strong 
fluctuations in this quantity (olq k y^l jl k 0.4) and consequently in the two-particle transfer cross 
section (cr ~ d^, i.e. fluctuation in a of the order of 100 %). 

In keeping with the analogy presented in Fig. [2| in the case of electromagnetic transition 
between members of a quadrupole rotational band one expects in heavy nuclei fluctuations of the 
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order of ( V250/250)^, i.e. less than 1%. Within this context the average value of the absolute 



experimental cross section reported in Table IV is 1551 juh, while the average difference between 
experimental and predicted values is 81 yub. Thus the discrepancies between theory and experiment 
are bound in the interval < \crexp(i f) - cfthii f)\/o-exp{i f) < 0.09, the average 
discrepancy being 5%. 



In Fig. 10 the excited, pairing rotational band associated with the average value of the 0"^ 
pairing vibrational states with energy < 3 MeV, is displayed together with the best parabolic fit. 
Also given is the relative (p, t) integrated cross section normalized with respect to the gs gs 
transitions, a value which is in all cases < 8%, in overall agreement with the single j-shell estimate 
(see ref. [lOJ App. H), given in the inset to the figure. The result testifies to the weak cross talk 
between pairing rotational bands and thus of the robust off-diagonal, long range order coherence 
of these modes. 



B. Pairing vibrational band in closed shell nuclei 



In Fig. 11 we display the expected pairing vibrational spectrum (harmonic approximation, 
see refs. [|T0l - [T2l and refs. therein) associated with the closed shell exotic nucleus ^^^Sn [|2l|3l, 
up to two-phonon states. Within this approximation, the one-phonon states are the pair addition 
\a) = IgsC^'^Sn)) and pair removal \r) = \gsC^°Sn)) modes. The two-phonon 0"^ (|;?i;(^^^Sn)) = 
\r) (g) \a) = |0+("^Sn); 6.5 MeV)) pairing vibrational {{2p - 2/z)-like) state of "^Sn, is predicted at 
an excitation energy of 6.5 MeV (see Fig. |3]). The absolute two-particle transfer differential cross 
sections associated with \a) and |r), namely 

^^^Snit,py^^Snigs), (£cm = 20MeV), (38) 
''^Sn{p,ty''Sn{gs), (£cm = 26MeV), (39) 

are reported in the insets. Using detailed balance the reactions above 

'^°Sn(?,;?)'^2Sn(0+;6.5MeV), (E cm = 20MeY), (40) 
^^^Sn(p,0"^Sn(0+;6.5MeV), (£cm = 26MeV), (41) 



are, within the harmonic approximation, equivalent to ( [38] ) and ( [39] ), exept for the relative flux 
which is determined by the ratio kf/ki. 
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Similar calculations to the ones discussed above have been carried out for the closed shell 
nucleus '°°Sn, the results being collected in Fig. 12 In this case the two-phonon 0^, pairing vi- 



brational mode of Sn is expected, again within the harmonic approximation, at an excitation 



energy of 7.1 MeV. As it emerges from Figs. 1 1 and 12 and at variance with the pairing rotational 



scheme, the two-particle transfer cross section associated with the excited pairing vibrational state 
is of the same order of magnitude than that connecting the ground states. 

Within this context, it could be intriguing to check whether the reaction 
^°^Sn(p,t)^''^Sn populates a 0^ state at an excitation energy of the order of 7 MeV. This state, 
can be written within the (pairing vibration) harmonic approximation, as |^*''*Sn(0^; 7.1MeV)) = 
\a) (g) \a) (g) \a) \r). Namely a four-phonon pairing vibrational state, where 

\a) = \gse'''-Sn)) ; |r) = Ig^^Sn)). (42) 

In other words, the reaction ^*'^Sn(p,t)^*''^Sn (0"^, 7.1 MeV) is, within the harmonic picture of pairing 
vibrations, equivalent to the reaction '*'°Sn(p,t)^^Sn(gi'). It is of notice that a three -phonon pairing 
vibrational states has been observed (53) in the reaction ^"'^Pb(f, p)^^^Fh at an excitation energy of 
about 6 MeV, with Q-values and absolute differential cross sections compatible with the excitation 
of the ^"^Pb pair addition mode, namely ^°^Pb(t,p))^'"Pb(g5'). While in this case deviations from 
the harmonic prediction are modest (essentially , most of them arising from the presence of the 
valence orbital pi/2 lying just below 6f(^°^Pb) 1541), in the case of pairing vibrations based on 
'""Sn, anharmonicities are expected to be much stronger. This is in keeping with the fact that 
N = Z nuclei display, as a rule, coexistence phenomena. That is, a strong competition between 
spherical and deformed 0^ states (cf. e.g. [|55l - [571l andrefs. therein, see also 11581 ). 



V. CONCLUSIONS 

The microscopic nuclear structure (BCS) description of pairing rotational bands, together with 
the second order DWBA description of two-nucleon transfer reactions which include successive, 
simultaneous and non-orthogonality channels provide, arguably without free parameters, an over- 
all account of Cooper pair transfer to superfluid nuclei. Inarguably, theory not only reproduces all 
reported '*"^^Sn(p, t)'^Sn(gs) absolute cross section data within experimental errors, but it does so 
with uncertainties below the 10% level. 

The study of the pairing vibrational scheme around the starting and end points of the pairing 
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Tables 



rotational spectrum promises to provide new insight on pairing fluctuations and their anharmonic- 
ities, in situations of large neutron excess and ofN~Z, i.e. around closed shell system ^^^Sn and 
of the, likely deformation coexistent '""Sn, respectively. 
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TABLE I: Results for the valence shell lying closer to Fermi energy, of a BCS calculation, for '^"Sn. The 
mean field used corresponds to a Saxon-Woods potential, in a 15 fm spherical box (continuum 
discretization). The pairing coupling constant used, G = 0.05 MeV leads to the value of pairing gap 
obtained from the three point formula (A = 1.47 MeV), summing the contributions of states from the I51/2 

at -39 MeV to -1-6O MeV. The resulting Fermi energy computed by solving the BCS number equation is 
A - -6.9 MeV. The quantity a'^ = i:^>o (= Zm>o U'jV'j = ^jU'jV'j) « 6.1 for the valence-shell 
space and a'^ w 29.4 for the whole single-particle space used in the calculation. 
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TABLE IF Results for the valence shell lying closer to Fermi energy, of a BCS calculation, for '^"Sn. The 
mean field used corresponds to a Saxon-Woods potential, in a 15 fm spherical box (continuum 
discretization). The pairing coupling constant used, G = 0.18 MeV leads to the value of pairing gap 
obtained from the three point formula (A - \ A1 MeV), summing the contributions of states from the I51/2 

at -39 MeV to -i-O MeV. The resulting Fermi energy computed by solving the BCS number equation is 
A - -6.72 MeV. The quantity ofg w 6.1 (see caption to Table|lj) for the valence-shell space and a'^ « 8.2 for 

the whole single-particle space. 
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TABLE III: RPA wavefunctions of pair addition and removal mode of Sn (above) and Sn (below). 
Single particle energies have been taken from experimental values referenced in the National Nuclear Data 
Center. The energy of the lowest pairing addition and removal phonons in ^^^Sn are respectively 
W{A + 2) - 3.45 MeV with G(A + 2) = 0.131 MeV and W(A - 2) = 3.06 MeV with G{A - 2) = 0.157 
MeV, the associated particle -pair vibration coupling strength (see Eq.([2T])) being Kadd - A(A + 2) - 1.08 
MeV and A^m = A(A - 2) = 1.60 MeV respectively. The minimum of the dispersion relation, and thus the 
Fermi energy, are equal to A = -4.75 MeV (see Fig. [3]^a)). The pair degeneracy of the single-particle space 

associated with ^^^Sn and ^"'^Sn is Q = 31 and £1 = 21 respectively. The energy of the lowest pairing 
addition and removal phonons in ^^Sn are respectively W{A -i- 2) - 5.13 MeV with G{A + 2) = 0.290 MeV 
and W{A -2) = 1.96 MeV with G{A -2) = 0.380 MeV, the Fermi energy being A = -14.5 MeV (see Fig. 
[3|b)). The associated A values being Aadd = 2.12 MeV and A^^ = 5.32 MeV. The binding energy of ^^Sn 
was assumed to be B(^*^Sn) = 794.24 MeV, from the polinomial (4th grade) fit of the binding energies of 

tin isotopic chain. 
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TABLE rV: Absolute cross section associated with the '*"^^Sn(;7, t)'^Sn(gs) cross sections (i.e. between the 
members of the Sn-ground state pairing rotational band) calculated as described in the text, in comparison 

with the experimental findings, 
jub; the number in parenthesis corresponds to the statistical errors; the numbers in square brackets 
provide the angular range of integration of the absolute two-particle differential cross sections, 
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FIG. 1: (a) The occupancy (Ny) ~ of the single-particle states of the unperturbed system, in which the 

individual pairspins are aligned along the gauge (z-)axis (see Fig. A. 1). This non-correlated system 
(a^ = 0), displays zero pairspin alignment {PS = 0), that is {S x) = 0. (b) The superconducting (nucleon 



superfluid) ground state displays Off- Diagonal Long Range Order (ODLRO, see Eq. ( A32 1) and a finite 
value of the total pairspin (PS 0; = Yjv>o U'yVy), i.e. {sx{v)} + 0, can be viewed as a one-dimensional 



domain wall (see also Fig. |A.3| of Appendix [A]). The quantities in the inset, are the amplitudes with which 
z- and x-components of the pairspin mix, leading to a privileged orientation in gauge space perpendicular 

to the gauge (z-)axis (see Fig. A. 1 (b)). 
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FIG. 2: Schematic representation of the nuclear structure consequences of spontaneous symmetry breaking 
of rotational and of gauge invariance (see also Table XI of ref. |[T2l \ 
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9.78 MeV 




M99/2) 



8.05 

2e(4/2) 



15.65 MeV 



FIG. 3: Dispersion relation (see ifTTll : see also ifTOl Eq.(5.50)) associated with the pairing vibration of 
^^^Sn(a) and of '°°Sn(b). (a) The part of the curve to the left of the minimum corresponds to the pair 

removal mode (|gi'('^*'Sn) >) while that to the right is associated with the pair addition mode {\gs{^^^Sn)}). 

The energy of the modes W{A ± 2) are measured from the minimum of the dispersion relation, its values 
being explicitly indicated (see also caption to Table 2). The coupling constants G(A ± 2) used in the 
calculations are given in the caption of Table 2, where the X and Y amplitudes of the corresponding 
wavefunctions are displayed. As seen from the insets, the pairing vibrational modes blur the sharp 

distinction between occupied and empty states. In fact, the pair addition mode can be excited not only by 
transferring two neutrons to levels lying above the Fermi energy, a process proportional to the 
Xaddik) = ^ '^-w^J'f'' amplitude (inset to the right), but also to states lying below the Fermi energy, a 
process proportional to Yaddii) = *" '^Iw^"'''' ■ The associated values of the particle-pairing vibration 

coupUng strength are Aadd = A(yS - +2) = 1.08 MeV and Arem - A{fi - -2) - 1.6 MeV. (b) The same as 
above, but for the closed shell system ^"^Sn. In this case Aadd = 2.72 MeV and A„m = 5.32 MeV. 
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FIG. 4: (a) Value of the U'V products associated with ^^"Sn (see Table 



II I. (b) Value of the XY products 



associated with the pair addition and removal modes of ^^^Sn (see Table 



nil 
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FIG. 5: Schematic representation, in terms of pairspin alignment, of the dynamical ODLRO induced in the 
ground state of the doubly magic nucleus ^^^Sn, by the zero-point fluctuations (ground state correlations) 

associated with pair addition and pair removal of the closed shell system (for comparison with the 
superfluid system '^°Sn see Fig. A.3 of Appendix [Xjl. The pairspin states are defined in Eq. (20ai (fy < ep) 



and in Eq. (20b I (ey > ep), the X and Y amplitudes corresponding to Eq. (21 1 (see Table 3). The mean 
square root value of the angle of the pairspin measured from the gauge (z-)axis is fixed by the relation 
cos^v - {sz{v))/\s\, which for occupied states is cosOy = (l/2)/(3/4)'^^, i.e. Gy - 54.5", the average angle 



of the precession cone of pairspins centered around the z-axis. Pairspin states (20ai and (20b i, have been 



calculated making use of the X, Y values reported in Table III and of 6y similar to the above one, but which 
for simplicity was chosen equal to be 45°. In the orientation of the pairspin reported in the figure, this fixed 
angle has been subtracted. In other words, the z-axis has been rotated by 45° into the z'-axis. Within the 
present scenario, it is expected that one-particle transfer reactions on '^^Sn may e.g. excite the 
® S^O^'^^^) ~ st^te of ^^^^Sn (see inset) with a weak, but likely observable cross section. 
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FIG. 6: (a) Radial function p(r) (hard core 0.45 fm) entering the triton wavefunction. (b) Radial function 

p(r) entering the deuteron wavefunction. 




FIG. 7: Predicted absolute differential '^'*'^Sn{p, t)'^Sn(gs) cross sections for bombarding energies 21 MeV 
< Ep < 26 MeV, and Ep = 40 MeV in comparison with the experimental data (see Il39ti44]l and ||45]| 

respectively). 
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FIG. 8: (Color online) Spatial structure of a two-neutron Cooper pair of ^^"Sn (see (|B16|),(|B17|) and 



(B23 1). The modulus squared wavefunction |^o(n? ^2)P = K6|n, ^2)^ (see Tables [l] and |n]), multiplied by 



ISn^r^r^ and normalized to unity, is displayed as a function of the cartesian coordinates xi = rjcosOu and 
X2 = r2sin0i2 of particle 2, for a fixed value of ri = xi = 5 fm (black dot) of particle 1, close to the surface 
of the nucleus (red circle). The numerical percentages correspond to the two-nucleon integrated density in 
a spherical box of radius 4 fm centered at the coordinates of the fixed particle. 
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FIG. 9: (Color online) Pairing rotational band along the tin isotopes. The lines represent the energies 
calculated according to the expression ^B = B(^°+^SnA,) - 8.124// + 46.33 1 10], subtracting the 
contribution of the single nucleon addition to the nuclear binding energy obtained by a linear fitting of the 
binding energies of the whole Sn-chain. The estimate of f?- Ill was obtained using the single j-shell 
model (see e.g. lITOll App. H).The numbers given on the abscissa are the absolute values of the 
experimental gs gs (in units of jub; see Table 4). 
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FIG. 10: (Color online) The weighted average energies (f'e.vc = Tji EiCTi/ cr,) of the excited 0^ states 
below 3 MeV in the Sn isotopic chain are shown on top of the pairing rotational band, already displayed in 
Fig. [9] Also indicated is the percentage of cross section for two-neutron transfer to excited states, 
normalized to the cross sections populating the ground states. The estimate of the ratio of cross sections 
displayed on top of the figure was obtained making use of the single j-shell model (see e.g. [J 0,1 App. H). 
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FIG. 11: Pairing vibrational scheme around ^^^Sn and calculated absolute reaction cross sections 
associated with the pairing addition and removal modes. The two-nucleon transfer spectropic amplitudes, 



used in the calculation are collected in Table III The optical model potential for the three channels, 
(namely A + t ^ (A + I) + d ^ (A + 2) + p or reversed) were taken from refs. 11411 and Il48l . 
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FIG. 12: Pairing vibrational scheme around ^""Sn and calculated absolute reaction cross sections 
associated with the pairing addition and removal modes. Concerning the spectroscopic amplitudes and 
optical parameters used in the calculations see caption to FigjTTj 
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Appendix A: Pair spin and domain wall 



In what follows we will discuss the mean field properties of the Hamiltonian 



where 



with 



and 



Hp-AN = H[p + Vp, 



v>0 



ey = Sy - A 



Ny = a^Gy + a^ay 



The pairing interaction is defined as 



V„ = -GP^P, 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



where 



and 



v>0 

p; = at a;. 



(A6) 



(A7) 



The single-particle Hamiltonian H'^^ is invariant under time reversal operations. As a conse- 
quence, orbitals are twofold degenerate, the corresponding states being denoted |v) and |v). 
Because the operators A^^, P^ and Py satisfy the commutation relations (cf. Appendix ICl) 



[P:,Py]=Ny-l 
[Ny-hPt]=2Pt , 



(A8) 
(A9) 



and 

[Ny-l,Py] = -2Py, ( A 1 0) 

one can define the x, y and z components of the pairspin operator s(v) according to the relations. 



1 

1 



' ^,(v) = ^(/':-/'v) = ^ 



-i 



(All) 
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(C) 




2 Z^i/>0 



FIG. A. 1 : (a) In the presence of the pairing interaction, the pairspin vector < ^ > acquires a component in 
the X, y plane. It is then possible to define an intrinsic system 7C', which is obtained rotating the laboratory 
system 7C by the gauge angle 2<p around the z- (gauge) axis (positive angles correspond to counter 
clockwise rotations). The contributions from the individual pairspins (all lying in the (z, x') plane) are 
schematically shown. The perpendicular component dominates for states close to the Fermi energy (e « 0), 

while the pairspins associated with states far from ep are aligned along the z-axis. (b) Contribution of a 
pairspin associated with single-particle states of energy ty. The pairspin vector makes an angle 6y with the 
gauge z-axis such that sin^y = A/Ey and cosOy = -ey/Ey. (c) The total pairspin vector is the sum of many 
individual contributions. The value of its projection on the z-axis is equal to (N - Q.)/2, where N is the 
number of particles and Q is the total pair degeneracy of the single -particle subspace considered to 



describe the system (see Eq. (B18 1). 
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and 



1 
-1 



(A12) 



In fact, using the commutation relations (A8)-(A10) stated above one obtains (cf. Appendix 



0, 



[sy{v), 5,(v)] = is,^(v) , 



(A13) 
(A14) 
(A15) 



Of notice that these s = (s^, Sy, s,) operators although acting in an abstract, gauge space, are 
as real as the standard spin of electrons and nucleons. The z-component of pairspin pointing up 
means "occupied" two-fold degenerate orbitals, pairspin pointing down means "empty", while a 
pairspin pointing sidewise implies a certain phased linear combination of up and down (see Fig. 
[T]as well as Fig. A.l I. In keeping with this scenario, the BCS ground state displays a gradual 
rotation, like a domain wall, of the pairspin vectors across the Fermi surface. 

The eigenvectors of s^iv) in pairspin space are 



|l)v = 



1 



(A 16a) 



and 



I2)v 







1 



= |0)v. 



(A 16b) 



To better clarify the meaning of state |1 >y and |2 >y, let us assume to be working with a set of 
two-fold degenerate states vi, V2..., each pair of levels connected by time reversal, e.g. (vi, Vi). In 
the uncorrelated case (G = 0), |1) and |2) can be viewed as fully occupied or fully empty states 



(see Fig. A. 2). That is, a two-particle (filled) and a two-hole state (empty) respectively (see. Eqs. 
(A. 17) and (A. 18) below; see also (A. 22) and (A. 23)). The same argumentation can be applied 
to each pair of (m, -m) states connected by time reversal, of a general set of (2j -I- 1) degenerate 
single -particle states. In other words, the system under consideration displays its pair addition and 
pair removal modes, at the level of individual pairs of time reversal states (v, v), building blocks 
of which Cooper pairs are built. Within this context, it is of notice that Cooper's model works 
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equally well if one thinks of it in terms of a correlated two-hole state in the Fermi sea, BCS being 
an extension and, in a way, a natural melting of the two views, as required by quantum mechanics 
(zero point fluctuations (ZPF) which, within the present context can be interpreted in terms of 
ground state correlations (gsc)). The quantal nature of these correlations is further evidenced by 
the fact that the diff"erent components enter the correlated Cooper pair in terms of probability 



amplitudes (see in particular (B 17 )). 

In the basis (A. 16), the operators s^, Sy and have the same matrix representation as the Pauli 
matrices except for a factor 1/2, that is, a = 2s. The action of and N on the states (A. 16) is 
given by 



^z(v) 



( \ 
1 







, 1 , 



1 







1 



2 


2 




1 

-1 

1 
-1 



1 







1 , 



1/2 



1 



-1/2 







, 1 , 







1 



1 



0. 



(A17) 



Inverting the relations (All), the pair operators and P can be identified with the raising and 
lowering operator in pairspace. In fact. 



= J](SX(V) + iSyiv)) =S, + iSy =S^, 



(A18) 



v>0 



and 



P = ^(5v(v) - isyiv)) = Sj,-iSy = S^. 



(A19) 



y>0 



In this space, i.e. the space subtended by the states |l)v and |2)y, the operators S +(v) and 5_(v) are 
represented by the matrices 



«/|5+(v)|j)) = 



<iip:ii> <iim2> 

<2|P+|1> <2\P:\2> 



\ 


( 









1 







0. 



(A20) 
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«/|5_(y)|7)) = 



< l\Py\l > < l\Py\2 > 
<2\Py\l > < 2\Py\2 > 





( 












) 


1 


0. 



(A21) 



while 



(OW(y)-llj)) 
Consequently 



while 







, 1 , 







, 1 , 



1 




( 

< yy |(A^y - 


Y)\vv > 


< VV\{Ny - 






( \ 

1 


^ < 0|(A^, - 


\)\vv > 


< 0\(Ny - 


1)10 > ^ 






\(\ ( 




( \ 









(A22) 







, 1 , 



, 1 , 



1 



0, 



(A23) 



= 



1 





1 



^ ( \ 
1 



( \ 




, 1 , 



(A24) 



The total pairspin in the z-direction is closely related to the number operator 



(A25) 



v>0 



where A^^ = 2v>()^v, and Q. is the total number of two-fold degenerated single-particle orbitals, 
associated with the single-particle space considered. Within this context, see Eq. ( |B18| ). 
The two terms of the pairing Hamiltonian can now be rewritten, as 



//;^ = J]6,A^, = J]6,(l+2^,(y)), 



(A26) 



v>0 



v>0 



and 



'z' = - Z ^^'r 



V2 



v\ ,V2>0 



-G 



Z Sxiyx) Z ^x(V2) + Z ^y(Vl) ^ Syiyi) 



yi>0 



yT>0 



v,>0 



vt>0 



v,>0 



(A27) 



vi,y2>0 vi>0 

where sj^) = Sx{v)i + Sy{v)j, i and j being unit vectors along the x- and y- directions. The 
last term is the contribution of the pairing interaction to the single-particle mean field Hsp. Al- 
though it can easily be incorporated in this term, it is customary to neglect it, in keeping with the 
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FIG. A.2: Schematic representation of the eigenvectors of ^^.(v) in pairspin space: states |l)v and \2)y can 
be viewed as a pair addition and a pair removal mode, at the level of individual pairs of time reversal states 

(v,v). 



schematic nature of V,,, tailored to act on the pair space. On the other hand, this contribution is 
to be considered when comparing the solution of schematic models with exact solutions (see e.g. 

m,m). 

The interaction Vp is a spin-spin coupling, which seeks to align the transverse components of 
the pairspins. Was it not for H^p, which depends on 5 j, all pairspins would line up in the same 
direction (strong coupling limit), perpendicular to the z-axis, in keeping with the fact that Vp is a 
function of only the S ^ and S y pairspin operator. In the opposite limit, that is, for Vp = 0, pairspin 
alignment is zero, that is, there is no component of the pairspin in the (x, y) plane perpendicular 
to the z-axis (see Fig. 1(a)). From this figure it is clear that H^p counteracts pair spin alignment. 
Indeed H^p plays, in the nucleus, the role of a magnetic field in a solid, which tends to align the 
spins in the z-direction, or opposite to it, with a strength (as measured by Sy ) that increases in 
absolute value as a function of the energy of the twofold degenerate levels (v, v) away from the 
Fermi energy sp = A. Thus, for sufficiently large values of \Sy - 6f\, of the order of 2 A (2-3 
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MeV), the single-particle energy as measured by Hs,, dominates (see Fig. 1 as well as Fig. lA.lTa) 



and A. 3 ). This is also the reason why calculations of pairing correlations in nuclei, which depend 
on small contributions arizing from orbitals distant from the Fermi energy are to be handled with 
care, in particular when comparing the calculated results with the experimental findings. Within 
this context see Sect 



IV A 



discussion connected with agCn, r2). Close to the Fermi energy, Vp is 



the overriding effect, and a large transverse pairspin is expected, as testified by the coherent sums 



in Eq. ( A27 ). While any single pairspin in these sums displays large quantal fluctuations, the total 
pairspin has a well defined magnitude and orientation, since the fluctuations of the constituent 
pairspins add quadratically. 

The pairing Hamiltonian can be diagonalized in the mean field approximation, substituting one 



of the sums in Eq. (A27) with its average value in the mean field ground state. Let us assume 
that the average value of the pairspin polarization vector K5*x)l = 'K'{BCS\S ^_\BCS)'K' (see. Eq. 



A51 below) is non-zero, and choose a definite orientation for it in the x - y plane (denoted x'), 
subtending an angle 20 with the x- axis (definition of the body-fixed, intrinsic frame, see Fig. 
A.l). The pair interaction is then replaced by a mean field, namely the pair field of strength 
A' = GK5^x)l- One can then write. 



Up = -2G\{^^)\ J] [^,(v)cos20 + Sy{v)sm2((>] 

v>0 

= -2A' ^ [5';f(v)cos20 -I- 5y(v)sin20j . 



(A28) 



v>0 



Making use of the relation given in ( Al 1 1 one finds, 

A'(^.v(v)cos2<^ + 5,(y)sin20) = ^ (pjA'^-^'^ + PA'e^"^) = ^(PjA + P^A*), 
where we have introduced A = e~^"^A'. Consequently, 

Up = -A' Y,(Py +Py) = - J](Pl^ + ^vA*), 



(A29) 



(A30) 



v>0 



v>0 



P'^, P' denoting the operators in the body-fixed frame of axis (x', y') (cf. Fig. A.l a) and Appendix 
B|) in which by definition (f) = 0, that is {Sy) = and (S'J = {S'^) = a[y 
The total Hamiltonian then becomes a sum over individual pairspins. 



(Hp)MF - ^ hy. 



(A31) 



y>0 
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where 





( \ 






1 




/Zy = 6y + 6y 




-A' 

V 



COS 20 - /sin20 
cos20 + hrnlcf) 



(A32) 



It is of notice that thinking in terms of the independent (quasiparticle) densities the first term is 
connected with the normal and the second with the so called abnormal density respectively, leading 
to ODLRO. 

Within this context, in the mean field associated with the (diagonal) first term of the above 
equation (see also Eq. (A. 2)), the particles which move independently of each other, are nucleons. 
In order that this can happen, all the nucleons must participate in a highly coherent ballet following 
a refined coreography, in such a way that each dancer moves as if he was alone in the scene, being 
fenced-in through a dancers-like wall, only when approaching the edges of the scene. As it has 
been stated in the literature ll6T1l . it is a "rather unfortunate perversity" which views independent 
particle motion as antithetic to nuclear collective motion. 



In the case of the mean field described by the second term of Eq. (A32), the entities which 
play the role of nucleons in the case above, are now pairspins (pair addition and pair removal (v, y 
modes)). The only circumstance in which pairspins feel the pushings and puUings of the other 
pairspins, is when they try to adopt a different orientation but that defined by 5x (i.e. .x:' -direction, 
see Fig. A.l), being forced to align back by a domain wall. In the present case, the ballet is not 
performed by single dancers in a scene, but by couples in a crowded dancing hall. In spite of 
such a less cultured setup, the coreography is even more refined than previously described. This is 
because the partners of each dancing couple can, not only when close to each other, but also when 
finding themselves at opposite extremes of the dancing hall (coherence length), follow the other 
partners moves without missing a single step. Such a coreography of strongly overlapping pairs 
translates, in the gauge (pairspin) space, into the definition of a privileged orientation. 

This can be better seen by expressing the diagonalization condition (|4]), namely 



const 



y>0 



v>0 



in terms of the quasispin operators ( Al 1 1 ( A12 1 and 



*^'(v) = -^(A^v-l). 



That is. 



A'^ 

{eyls-Xv) - A'25,(v)] + _^(6y -Ey) + — = -Yj EylsAv) + const. 



(A33a) 



(A33b) 



(A33c) 



v>0 



v>0 



v>0 
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In a similar way in which 5j.(v) is diagonal in the independent (pair) particle situation, Sf{v) is di- 
agonal in the quasiparticle representation, with eigeinvalues 1 /2 when acting on the corresponding 
occupied states (\HF) and \BCS) states respectively), and -1/2 when acting on the correspond- 
ing unoccupied states (levels above Sf and two quasiparticle states respectively). Equating the 
(quasispin) operator terms and the c-number terms one obtains 

A'2 



const 



= _^(6v - £v) + 



(A34a) 



v>0 



and 



€y A' 

— 5z(v) + — *.v(v) = S,'(V). 



(A34b) 

The above equation determines the angle 6y in the z - x plane, which leads to independent 
quasispin motion. In the strong coupling limit (|A'| >> \€y - A\), "magnetization" is total, all 
pairspins (pair addition and removal (v, v) pairs) pointing along the x- direction. 

As already stated in connection with Eq. (|4]) of the text, the c-number const is equal to the 
ground state energy, also known as the U term of (Hp)MF (see e.g. [fTOl Eq. (G.ll) of App. G; see 
also below, subsection on ground state energy and pairing correlation energy). 

In keeping with the fact that the quasispin states are normalized, and that one has chosen a 



representation in which 5„(v) = 0, the quasispin prefactors of Eq. ( A34b I must fulfill the relation. 



- 1 , 



(A35) 



It is of notice that the eigenvalue equation ( |A32[ ), 



fv - 



\ 


( \ 

Vy 


= {ey 


-E) 


( \ 

Vy 




Uy 









(A36) 



has the solutions +Ey, where (see also (A35 1) 



Ey= ^el + A'2 



(A37) 



The positive sign corresponds to the lowest total pairspin energy (6v - Ey), while the negative 
sign corresponds to the excited states (2-quasiparticle states, cf. subsection below on excited 
states). The relation between the Uy and the Vy components of the eigenvector can be deduced 



from equation (A36), 



-eyVy + A'e-^''^Uy = EyVy 
+ eyUy = EyUy. 



(A38) 
(A39) 
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Eq.( A38) leads to {Ey + ey)Vy = A'e ^"^i7y, implying a phase difference -20 between the Vy and 



the Uy occupation amplitudes, This allows one to write the following relations 



Uy = uyt 

Vy = V:e-'f 



(A40) 
(A41) 



with U'y and being the moduh of Uy and Vy respectively. These quantities can be calculated 



from the square modulus of Eq. ( A38 1 



(ey + EyfWyf = /^^\Uy\\ 



(A42) 



Making use of the normalization relation \ Vy\^ = 1 - \ Uy\'^, one obtains 



(A43) 



Defining the angle 9y according to 



cos^v = 



sin^v = 



A' 



(A44) 



one can rewrite the quasiparticle amplitudes as 



U'y = sm(ey/2) , K = cos(0v/2), 



(A45) 



where the angle 9y represents the angle between the direction of s(v) (z'-axis in Fig. A.l ^b)) and 
the z-axis. The average value of s^iv), Sy(v) calculated with these eigenfunctions, generalizations 



of the eigenvectors introduced before (see Eqs. (A16a) and (A16b)), are given by 



-iUyVl + VyUl) 



sin^ cos^ cos20 = U'yV'y cos20 = ^— cos20, 
2 2 2£v 



{UyVl-VlVy) 



sin— cos— sin20 = U'^'y sin20 = sin20, 
2 2 2£'y 



(A46) 



(A47) 
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FIG. A.3: Pairspin distribution associated with the seven valence single-particle orbitals lying around the 
Fermi energy of ^^"Sn. The calculations were carried out making use of the results displayed in Table [ill i'^ 
particular the values l]'yV[„ following the prescription discussed in the caption to Fig. |5j 



and 



= -{VyV:-uiu;) 



2E, 



(A48) 



The dependence of (^v(^)) and {Sy{v)) on the the angle (p is consistent with the ansatz (A28). 



Furthermore, making use of Eq. (A47), one can determine the modulus of 5x selfconsistently 



Y,^Uv))\ = \{s^)\ = h^y = Z ^^^^ = < 

v>0 v>0 y>0 



(A49) 



a relation which is closely connected with the BCS gap equation. In Fig. A.3 the pairspin distri- 
bution associated with the valence orbitals of '^"Sn, calculated making use of the results collected 
in Table [n] and the scheme discussed in the caption to Fig. |5j is displayed. 

The wavefunction describing the ground state of the system is the product of all pairspins, 
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each of which is a linear combination of the pair addition and removal (v, v) modes with the 
corresponding weights Vy and Uy respectively, amplitudes which define their alignment in gauge 
space. That is, 



f](f/,i2)v+Kii)v)=n 



v>0 



v>0 



/ ( \ 



Uy 

1 , 



1 



Y[{Uy + Vyala^^\0) 



v>0 



= W [e'l'V, + e-''^V:alal) |0) = e'"^ f] ([/; + Ke'^'^a^^) |0). 



(A50) 



y>0 v>0 

Leaving out the overall phase one can write 



|5C5(0)k = n (^^ + Vy^'I'ata;) |0) = 

v>0 

Y] (U'y + Vya:d-:) 10) = |5C5(0 = 0)k<, 



(A51) 



v>0 



where 7C and TC' denote the laboratory and the body-fixed, intrinsic frame in which by definition 



has = (see Fig. B.l I. 

The creation operator in the intrinsic, body-fixed frame of reference is 



(A52) 



where ^ = e is the gauge operator inducing rotations in the two-dimensional gauge space. 
The quantities U[, and Vy are real. 

In the independent particle limit, that is. 



1 



v>0 



V>0,y<V;r 



det(ll,22...A^A^)|0), 



(A53) 



as expected. 

Summing up, the instability of the Fermi surface associated with transverse pairspin polariza- 
tion is associated with a many-body wavefunction, product of the individual pairspin states. 



|0>V = Uy\S,iv) = -1/2) + Vy\S,(v) = +1/2 >= Uy\2 >y +Vy\l >y. 



(A54) 



superposition of pairspin up and down and therefore non axially symmetric with respect to the 
gauge axis (z-axis). In other words, a linear combination of pair addition and substraction modes 



mixed at the level of individual pairs of time-reversal states (v, v) (see Fig. A.2). In the same way 
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as r^Fh(gs)) spends part of the time in the state \™Fh(gs)) and part in \^^^Fh(gs)), r"Sn(gs)) is a 
mixture of pair addition and removal states |1) and |2), respectively. The transverse polarization - 
which inherently breaks gauge symmetry - arises from such a superposition, a phenomenon which 
is produced by the action of the "external" mean pair field Up. 

The pairspin polarization may rotate collectively about the z-gauge axis. The azimuthal angle 
is therefore a dynamical variable (pairing rotations). The static pair field constitutes a deformation 
that defines an orientation. Through this deformation, the system spontaneously breaks away from 
axial symmetry, and the indeterminacy in the number of particles in a pair correlated state is an 
inherent feature of this symmetry breaking. The static deformation introduces a collective degree 
of freedom and gives the system the ability to rotate as a whole around the gauge axis. 

When the mean field solution leads to (5^^) = 0^ the intrinsic motion has axial symmetry 
and hence conserves particle number. In these systems, gauge invariance can also be broken 
dynamically, a phenomenon which gives rise to the pairing vibrational spectrum observed around 
closed shell nuclei, in terms of highly enhanced, single Cooper pair tunneling processes. 



Excited states 



The BCS equation in the body-fixed frame. 



6, -A' 
-A' -6, 



V 



= -E 



V 



(A55) 



leads to the eigenvalue equation 



(6y + E)i6y -E)- A'2 = 0, 



(A56) 



with eigenvalues Ey = yje^ + A'^, corresponding to the ground state previously considered (cf. 



Eq. ( A37 1) and -Ey, corresponding to excited states. The associated eigenvectors are 

Igs) = 



|exc) = 



-V 



(A57) 



The state |exc) is a two-quasiparticle state. It can be excited acting with UyU^ on \BCS)'k. In fact, 
using Uy = U'yQ'y " A ^ ouc finds (it is of notice that s[{v) = s^(v)) 



a^a; = U'y^P': - V'y^P'y + 2 U ^ ^iy) = 



U'V U'^ 



_y'2 -U'V 



(A58) 
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and 



a,a, = U'JP' - V'^P': + lU'V's'Xv) 



( \ 
U'V -V'^ 

V V V 



u:f -c/:k 



(A59) 



Making use of Eqs. (A23 A24), one finds that 



Qr>t|gs) = Q'>i 



v:. 



U'V u 



r2 



_y'^ -U'V 



U', 



-V 



|exc), 



(A60) 



and 



Analogously, one finds 
and 



u:. 



-V 



U'V U'^ 
_y'2 -U'V 



u:, 



-V 



= 0. 



a^aylgs) = 0, 



ayOTylexc) = Igs). 



We also remark that 



2is'y{v) = P'; - = a>+ - a,a,. 



(A61) 



(A62) 



(A63) 



(A64) 



One can then show that the operator 2isy{v) acting once on the ground state produces the excited 
state and, acting twice gives back the ground state, but for a sign change, that is. 



2/5'(v)|gs) = |exc) , 2z5'(y)|exc) = -|gs). 



(A65) 



Thus, two-quasiparticle excitation is equivalent to a rotation of the intrinsic system by an angle n 
around the y-axis. This is in keeping with the fact that e^'^" = 2iS y. 

Ground state energy and correlations 

The energy of the ground state of the total system at the mean field level is obtained as the sum 



of the energy of each pair spin, - Ey (cf. eq. (A36)), which is the eigenvalue diagonalizing the 



terms eyNy - (P^A + PA*) of the mean field hamiltonian, plus the constant term ^ (cf. eq. (|4|)), 



^ A'2 

Y,iey-Ey)+—. 



(A66) 



v>0 
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In the case of the non-interacting system (A 0) Ey ^ \ey\, so that the system's ground state 
energy reads 

= Z^^^ - i^^i^- ^^^^^ 

v>0 

Consequently, the contributions corresponding to 6y > (i.e. fy > A), cancel out, remaining only 
those of the states below the Fermi energy (ey < 0). Thus, 

El = 2 ^ ey, (A68) 

y>0;e,,<0 

in keeping with the fact that in the present case all the pairs below the Fermi energy fully occupied. 
An important quantity characterizing superfluid systems is the so called correlation energy, E^orr- 
It is defined as the difference between the energy of the interacting system and that of the non- 
interacting one, 

A'^ 

E,.orr = Eg,-El = Y,iey-E,) + — -2 ^ 6y. (A69) 

v>0 v>0;ey<0 

Collecting the contribution arising from levels displaying 6y > and 6y < 0, the above equation 
can be written as 

A'^ A'^ 

Ecorr = Yj i^y- Ey) + Yj ("^^ ~ + = Z ~ + Z d^vl - ^v) + — • (A70) 

v>0;e^>0 v>0;e,,<0 y>0;e„>0 v>0;e„<0 

Assuming a symmetric distribution of levels around the Fermi energy, the above expression be- 
comes 



A'2 

E I ] 

v>0;ev>0 ^ 



corr = 2 Y (ey-£y)+— . (A71) 



It proofs useful to express the ground state energy in the term of occupation probabilities and 
potential interaction energy, namely as (see term U eq. (G.l 1) of ref. [10] ) 

E,, = Y^eyV:'-^. (All) 

v>0 

Using the expression for V^^ 



y>0 \ 

which may be rewritten as 

A'2 



y>0 ^ ' ^ v>0 y>0 ^ ' 

^ ^ A'2 A'2 ^ A'2 A'2 



G 



Ey G 

y>0 y>0 ^ y>0 

A'2 



= Y(^y-Ey)+—, (A74) 



G ' 

y>0 
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which coincides with (|A66|). 



Appendix B: Pairing rotational band wavefunction 



X : laboratory frame) 




CK' : body fixed, 
intrinsic system) 



FIG. B.l: Schematic representation of the deformation in gauge space associated with a superfluid nucleus 

leading to pairspin alignment (see also Fig. 



A.li 



Let us start by defining the operator inducing a gauge transformation. In keeping with the fact 
that the generator of such a transformation is the particle number operator, we will be dealing with 
rotations in a two-dimensional space. 

The operator inducing rotation in pairspin space about the gauge z-axis is given e'^"''''^^^'^ The 
representation of this operator in pairspin space is 



e"'^ 
e''^ 



(Bl) 



This operator converts the state \(p = 0)y = 



v: 



, into the state rotated by an angle (cf. Eqs. 



(A40MA41)): 



e-"^ 
e'^ 



v: 



(B2) 
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In the following we shall use the operator ^y(0) = e which differs from e ^'^'"^^^ only by an 
overall phase e"^. Its representation in pairspin space is given by 



g-2i^ 

1 



(B3) 



and its action on 



>, 



is given by 



1 



v 
u: 



u:, 



(B4) 



producing the phase difference characterizing the rotated, |5C5(0) > state (cf. Eq. (A.51)). 
Calculating the average value of P in the (rotated) state Q{(p)\4> = 0)y, that is ao(v) = y{(p = 
0\g-\(l))Pg{(P)\(P = 0)y, one finds 





1 



U'V. 



(B5) 



We can also calculate the same average value by rotating the operator (Heisenberg representa- 
tion), rather than acting on the state. The rotated operator is given by Ql{(t))Py{(p = Q)Qv{<P)- The 
action of Ql{^)Py{(p = 0)Qy((f)) transforms the initial operator in the intrinsic frame Py((f) = 0), in 
which it has the average value a'^iv) = U'yV'y into the laboratory frame, in which its average value 
is aQ{v) = a'f^(v)e~^'^ . The inverse transformation, from the laboratory into the intrinsic frame, is 
effected by the operator {Ql{(p)Py{(l))Qy{(p)y^= Qy{4>)Py{<l))Q\4>). In what follows, we shall use 

P'y = Py{<P = 0), and Py = Py((f>). 

Similar considerations can be applied to the many-body wavefunctions. Let us consider a 
wavefunction 

= a{4...a]yj0 >, 

with a fixed number of particles. Let us now apply to the creation operator 



(B6) 



a'^ being referred to the intrinsic, body-fixed reference system (see Fig. B.l ). Thus 

¥^(<^) = e''''fa'^a'l..a'l\0 >= e'^^WK'icP = 0), 



(B7) 



(B8) 
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In keeping with the fact that 



O0 



d(p 



(B9) 

(BIO) 
(BU) 



This is equivalent to saying that 



( d \ d dil/ dil/ 

i.e. [0, A/^] = i. The BCS wavefunction can then be written, in the intrinsic body-fixed frame, as 



\Bcs((p = 0) ><K'- Y\ ^^1^ >~ n ^^^^^^ ^ ■ 

V v>0 

The corresponding normalized wavefunction is then 

\BCS((P = 0))<K' = Y\(U: + Vyja'^m = Y\(K + e-''Xalal)\0 > 

v>0 V 

( g^2(<^ -4i<l> ( ^ 

= l5cs(0)k = ([]f/;) i + ^J]cX4 + ^ Yj'^A^] +■■■■ 



(B13) 



v>0 



v>0 



Vv>0 



(B14) 



where Cy = Vy/U'y. Thus, 

\No) = J d(l)e'^'^\BCS{(l) = Q))'K' = J d(pe'^"^\BCS {(p))^^ 



= (n,>of/;) I ^^'•^''^(i + ... + 



-iNcj) 



(N/2)\ 



2 ^y^t^t 



\N/2 



^v>0 



+ ...)I0) 



.No/2 



YjCyulaU 10). 



(B15) 



. v>0 



It is of notice that the factor e''^"^ above is equivalent to the transformation coefficient e'^^ between 
p and q representations. Finally, the members of a pairing rotational band are described by the 
states, 

/ \A'o/2 

' (B16) 



2 Cyalal 10). 

v>0 / 

Now, from this relation it becomes clear that the Cooper pair wavefunction 

\6) = Y,Cyalal\0), 



(B17) 



v>0 
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is to be interpreted to be valid for values of Sy close to e^, otherwise one risks not to be able to 
normalize it {Uy for Sy « ep < for deeply bound occupied states). This is in keeping with 
the fact that pair condensation in general, and nuclear superfluidity in particular, are associated 
with a modification of the Fermi surface within a narrow band around it (sf + A). In other words, 
pairspin alignment as described by BCS implies that something unique takes place in the long 
wavelength limit of the spectrum, namely the appearance of a coherent state with almost (aside 
from the weak (Sy « Sf) dealignment introduced by H^p), perfect phase coherence. This state 
behaves essentially semiclassically, and its properties can hardly depend on the ultraviolet behavior 
of the system. In other words, Ecuwff can be set to include only the valence single-particle shells, 
adjusting G to reproduce the value of the pairing gap. Within this context see also the discussion 
in Sect. 4 in connection with Fig. [8] and Tables |l] and |llj 
In the case of a single j-shell (see e.g. [ITOl App. I), 



V = 




U= J\ 



N 

2Q' 



thus 



while 



U'V = 



2^^ I 2D.) 



U 



1 



N_ 
2S1 



N 



2Q-N 



(B18) 



(B19) 



(B20) 



For a number of particles considerably smaller than the full degeneracy of the single-particle sub- 
space in which nucleons can correlate, that is for << 2Q, one can write 



U'V = 



N 

2h 



1 



A^ 
4Q. 




and 



Consequently 



U' 



N 



2Q. 



(i-i) 



A^ 



A^ 



2Q. ^ 4Q 



A^ 

2h 



U'V 



10) « UyVyalallO), 



(B21) 



(B22) 



(B23) 



v>0 



in keeping with Eqs. ( 15b I and (28 1. 

In what follows we work out some relations which are useful to calculate expectation values in 
the |5C5) state. 
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Making use of (B7) (i.e. a'J = QunQ ^ = e "^al) one can write 



thus 



Making use of ( lb I and ^ one can write 



a , 



't 



0" = y). 



VX 0« = v), 



(B24) 



(B25) 



(B26) 



in keeping with the fact that U'^ and are real c-numbers, and thus U'^ = Up and V'^ = V^, and 
consequently Uy = Uy and Vy = Vy, as well as the fact that = -ay, a consequence of the anti- 
unitary character of the time reversal operator. Within this context it is of notice that the intrinsic 
property of a nucleon of being in a state with quantum numbers ( j, m) or ( j, -m) does not of course 
affect the gauge angle of rotation (20, see also Fig. |A. 1 1 defining the intrisinsic (body-fixed) frame 
of reference with respect to the laboratory system. Taking the hermitian conjugate of the second 



case of eq. (B26) one obtains 



(B27) 



Multiplying the first entry of Eq. ( B26| ) by Uy and (B27 1 by Vy one obtains, 



U'yal = uy: - V'yUX^ 



V'ya-y = V'yU'ya'y ^ Yy^d^ 



Summing these two expressions and making use of ([3]) one obtains 



a} = U'ya\ -I- Y[,ay 



which is equivalent to 



Uya\ + Vla-y. 



(B28a) 



(B28b) 



Similarly, multiplying the hermitian conjugate of the first entry of Eq. (B26 1 by V[ and the second 
entry by U'y and subtracting the resulting expressions leads to 



d^, = u'A. - v[ 



ay 



(B29a) 
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which is equivalent to 



at = Uyal - V*ay 



(B29b) 



One then obtains 



= Yj ^Wv = Z l^yy'^l - (K)'«vffv - UyV:(ala, + ala,) + f/,K) , (B30a) 



v>0 



v>0 



v>0 



v>0 



thus 



ao = {BCS\P\BCS) = J] = ^""'^ J] ^^^^ = ^""^^o 



v>0 



v>0 



2 f^vK 



Vv>0 



{BCS\P^\BCS)* 



Summing up 



and 



ao = ^ UlVy, 

v>0 



v>0 



leading to 



A = Goto = e'^'^Ga^ = e'^'^A'. 



(B30b) 



(B31) 



(B32) 



(B33) 



(B34) 



An alternative derivation of the above relations can be obtained by inserting the expressions of d 



't 



and a'J obtained from (|B28a|) and (|B29al) into 



P = J^ a,ay = e'^'* Yj = e^^'I'P', 



(B35) 



v>0 



v>0 



which leads to 



P = YjiUla, - y>I)(f/>v + V^al) 

v>0 

= [VyC^v^v - V'yaWy - Uy,{alay + ala,) + U'X] . 



(B36) 



v>0 
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Appendix C: Commutation relations 

Making use of the relations 

[AB, C] = A[B, C] + [A, C]B, (CI) 

[C, AB] = A [C, B] + [C, A]B, (C2) 

[AB, C] = A{B, C} - {A, QB, (C3) 

and of the definitions 

= a^a^ ; Py = a^ay, (C4) 

one can calculate the commutator 

[Py,Py] = [ayay,a'la^'\ = UyiOyay, a^~\ + [ayay,a^'\a^ (C5) 

= (-{ay, a^}ay) + (ay{ay, a^}) (C6) 

= -a^Uy + QyOi, = 1 - {a^Qy + Qy) = I - Ny, (C7) 



where 



Similarly, 



and 



Thus 



Consequently, 



Summing up. 



A^v = UyUy + a^tty. (C8) 

[d^ay,P^] = [ayay,a^a^] =a^[a^ay,a^] + [a^ay,a^]a^ = (C9) 

{uy, a*} = a^a^ (CIO) 

[a^Oy, P""] = [a^Uy, a^a^] =a^[a^ay, at] + [a'^Oy, a^]a^ (CI 1) 

= a>t{av,<}) = flX • (C12) 

[Ny,P:] = 2P:. (C13) 

[Ny,Py] = -2Py. (C14) 

[P:,/',]=A^,-1 (C15) 
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[A^,-i,p;] = 2p; (C16) 

and 

[N,-\,PA = -2P,. (C17) 
Making use of these relations and of the definitions 



^x(v) = -(P; + Pv) , syiv) = -{Pt - Py), (C18) 



and 



one obtains 



^z(v) = ^(A^v - 1), (C19) 



[^.v(v), Syiv)] = [^-(P: + Py), ^XP: - Py)] = ^(-[^v ' ^v] + [Py, P^) 

= hPv, PU = - A^v) = is,(y), (C20) 
2i 2i 



[Syiv), S.iv)] = ^XiPt, (Ny - 1)] - [Py, (Ny - 1)] = 'l-pPt + 2Py) 

= -^XPt + Py) = is.(v), (C21) 



[S,(V), S,(V)] = ^[{Ny - 1), Pt + Py] = ^((Ny " 1), P^] + [{Ny " 1), Py]) 



= -^{Py-P:) = +iSy(y). (C22) 



Appendix D: Generalized Rigidity in Gauge Space 

Generalized rigidity in gauge space implies that if one pushes, with the help of a field that 
changes the number of particles in two, one of the poles of a deformed system in gauge space 
(see e.g. Fig. B.ip , system which can be viewed as a wavepacket in particle number, the other 



pole reacts rigidly to the push, and the system starts rotating as a whole (pairing rotational band). 
Similarly, if it was already in rotation it changes its rotational frequency from to = A{N())ITi to 
oj' = A(No + 2)/fi, A being the Lagrange multiplier which, in BCS theory, is closely connected with 
the particle number equation. 

If the gauge space image is not sufliciently concrete to create a physical picture of the process, 
let us think of a quadrupole deformed nucleus whose intrinsic state is described in terms of the 
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Nilsson intrinsic state. Making use of a proton beam which acts upon one of the poles, the system 
reacts as a whole and starts rotating with a frequency associated with one of the allowed values of 
the angular momentum, the transfer quantum corresponding to an energy inversely proportional to 
the moment of inertia. 

One may argue that the reaction of the pole not acted upon by the external field is not in- 
stantaneous but takes place only after an interval of time, compatible with the propagation of 
information in the nuclear medium, has elapsed. This parlance is not even wrong, as a wavefunc- 



tion, in particular that describing the intrinsic ground state of a superfluid nucleus (see Eq. ( 12)), 
is not a matter function but a probability amplitude function [|62ll with perfect phase coherence 
throughout. Within the quadrupole deformed nucleus analogy, generalized rigidity implies that 
the nucleus reacts rigidly as a whole to the action of the proton field acting on a pole [|63l , even 
if the moment of inertia of the associated rotational band is that of superfluid nuclear matter, and 
thus considerably smaller than the rigid moment of inertia. 

The specific experiment to study the consequences (emergent properties) resulting from a spon- 
taneous breaking of symmetry (e.g. of rotational invariance) is a probe which itself violates the 
symmetry in question. Now, while most of the devices we find in a well equipped nuclear labora- 
tory violate rotational invariance - think for example of a proton beam line defining a privileged 
orientation in 3D-space - one does not find many which violate gauge invariance. In other words, 
while rulers and goniometers defy empty space isotropy and homogeneity, one does not usually 
walk around with instruments which do not have a fixed number of particles. 



This was the real importance of the Josephson effect (see Fig. D.l ), which provided a simple, 
and quantitative accurate answer to the question: how does one measures the gauge phase of 
a superconductor? The answer is, with the help of another superconductor displaying also an 
unknown but nonetheless well defined gauge phase. Establishing a weak coupling (oxide layer) 
so that electrons can tunnel, one at a time, between the junction. If the first system can be viewed 
as a rotor in gauge space, the second one can equally well be represented in this way. Biasing the 
junction with a constant potential difference, will lead to a two-rotor-coupled system (through pair 
transfer across the junction), rotating with frequencies which differ by ebMlfi = {Ai - A2)/fi. Such 
a system will display a resonant behavior (alternating current with frequency eAVt/h) provided 
the rotor, deformed system picture in gauge space, is applicable. The fact that this effect provides 
the most accurate measure of (e/fi) available, testifies to the validity of the deformed rotor picture 



in gauge space associated with the BCS wavefunction (12). 
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It is sobering that this is so, in keeping with the fact that it was Bardeen the most strenuous 
opponent to Josephson's ideas, even more so if one is reminded that the choice argument of such 
opposition was based on the fact that the pairing gap A = Gao vanishes at the junction. Now, the 
order parameter of BCS theory (as well as the justification of Cooper pair model) is otq = (P'^) = 
(P)*, namely the condensed pair field. Electrons may tunnel one at a time, without obliterating the 
validity of pair transfer, as the coherence length = fiVf/2A) is much larger than typical junction 
dimensions, diverging at the junction (kind of extreme ^ 'Li-halo-like phenomenon (see e.g. [|64l . 
Il65ll ') within the condensed matter framework). 

Within this scenario one can posit that, in the nuclear case, the equivalent of the Josephson junc- 
tion device or better , the setup of an experiment which can measure differences in gauge phases, 
allows for embodiments which are not possible within the field of condensed matter physics. This 
is related to the fact that fluctuations in finite many-body systems, in general and pairing vibra- 
tions in normal nuclei in particular, are not only quantitatively but also qualitatively stronger than 
in condensed matter. 



Consequently, not only a collision between two superfluid nuclei (see Fig. D.Lb)) can be 



viewed as a time dependent generalization of a Josephson junction (see. Fig. D.l 'a)). Also 
a normal- superfluid nuclear reaction can provide similar information, as a closed shell system 
displays a very collective pairing vibrational spectrum which can be viewed as large amplitude 



dynamical gauge symmetry violating mode (see Fig. D.l c)). 

Let us consider the ground state of a (light) closed shell system. Within the present discussion, 
it can be written as 

\^)) = a\aigs)) + fi\{a - 2){gs) ® (a + 2)igs)) (Dl) 

where \(a-2)(gs)) and \(a + 2)(gs)), are the pair removal and pair addition modes, while = 1 

is the normalization condition. In other words, the closed-shell system is part of the time in states 
with two more (see (I) inset Fig. D.l[ ), or two less correlated nucleons (see (II) inset Fig. D.l). The 
two-hole uncorrelated states which, arguably, ensure particle number correlation, in fact describes 
the ground state correlations (backwardsgoing amplitudes, see e.g. Table [III] Y-amplitudes as well 
as insets in Fig. |3]) associated with the pair addition mode which, dynamically, deforms the nucleus 
in gauge space defining a transient, privileged orientation. The same can be said concerning the 



two-particle uncorrelated system shown in the inset (II) of Fig. D.l c). 

Let us now return to the analogy with quadrupole rotational and vibrational bands. Going away 
from closed shell nuclei in medium heavy systems the energy of the first 2"^ state lowers in energy, 
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the corresponding period becoming longer, the associated amplitudes larger (see e.g. IfTDl Ch.7). 
Eventually, after the quantal phase transition has taken place, the rotation of the system as a whole 
can be viewed as a very low frequency quadrupole vibrational mode (for which the restoring 
force vanishes while inertia remains finite), which dynamically defines a privileged quadrupole 
deformation symmetry axis (and thus an associated set of Euler angles), orientation which after 
each period changes orientation, with a frequency inversely proportional to the inertia of the mode. 



In a similar way the state (Dl I defines dynamically, a privileged orientation in gauge space 
which specifically can probe the corresponding static quantity of a superfluid target nucleus A (see 



Fig. D.lc)), thatis. 



a+A^l a((a + 2) O (a - 2)(gs)) + A 



V N 



{a - 2){gs) + (A + 2){gs), 
(a + 2)(gs) + (A- 2)igs), 



(D2) 



where A (Xw cn\N) labels the ground state of a superfluid nucleus, e.g. of '^"Sn which can be 
viewed as a wavepacket in neutron number, the scattering state within curly brackets being a 
virtual set of states each displaying a dynamical or a static privledged orientation in gauge space 
and thus a gauge phase. A particular embodiment of such a reaction can be 



3Li6 + 50 ^1^70 



3Li4 + 5Q^Sn72, 



j'Lig + jQ^Sngs, 



(D3) 



in keeping with the fact that N = 6 corresponds to a (parity inversion) magic number, fhiigs)) and 
l^^hiigs)) being the pair removal and pair addition modes of ^Li (see [|66ll67l ). 



Appendix E: Kramers-Kronig dispersion relation 

The nuclear superfluid phase, and associated pairing rotational band behaves, because of its 
ODLRO essentially as a classical (coherent) state. This is the basic reason which is at the basis 
of the results displayed in Fig. [7] and Table IV that is the remarkable quantitative accuracy with 



which theory provide an overall account of the absolute value of the experimental findings. 

There is, however, a second reason for the higher accuracy with which one can, in principle, 
predict absolute two-nucleon transfer cross sections, as compared with one-nucleon transfer cross 
sections, or, simply, (although, arguably, less well defined see e.g. (511681 169||) absolute values of 
the single-particle spectroscopic factors. This is connected with the ambiguity and eventual lack 
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of consistency of empirically determined optical parameters. Typical of the first type of limita- 
tions, is connected with the fact that the depth of the real part of the optical potential can have 
difl'erent, commensurable values, all leading to the same phase shift (e.g. socalled Igo's ambigu- 
ity [TTOII ). Concerning the second point (consistency) , one is reminded of the fact that the real 
and the imaginary part of the optical potentials (U + iW) controlling off"- and on -the energy shell 
processes, are the real and imaginary part of the nuclear mass operator (sum of polarization and 
correlation contributions), referred to, also, as the nuclear dielectric function [[51 . Consequently U 
and W must fulfill the Kramers-Kronig dispersion relation (cf. refs. [fTTl [721 . see also [(H). This 
is a bidirectional mathematical relation, connecting the real and imaginary parts of any complex 
function that is analytic in the upper half of the complex plane. The Kramers-Kronig is a rather 
fundamental relation, in that it is strictly related to causality. Summing up, the real and imagi- 
nary parts of the optical potentials empirically determined from a global elastic scattering fitting, 
should respect the above mentioned dispersion relation. This condition is only marginally fulfilled 
in a number of cases, thus introducing uncertainties difficult to control. In keeping with the fact 
that, theoretically, U (r) results from the convolution of the nucleon-nucleon interaction with the 
nuclear density p(r), and that this quantity is best known around the nuclear surface, in a similar 
way in which W{r) receives important contributions from the interweaving of single-particle mo- 
tion and collective surface vibrations, one can posit that the lack of consistency between U and 
W are likely to be more serious for values of r larger and smaller than the nuclear radius (i.e. 
r > Rq and r < Rq). The phase coherence of the nuclear condensate implies that Cooper pair 
wavefunctions, closely related to the eff'ective two-particle transfer nuclear formfactors, are rather 
compact objects, essentially concentrated on the nuclear surface (see e.g. Fig. [8]). Consequently, 
the main contribution to the absolute two-particle diff"erential cross section arises from values of 
r X Rq (nuclear surface), region in which the optical potential is best known. On the other hand, 
the formf actor associated with single-particle transfer is a standard mean field wavefunction (in 
any case as far as the main peak of the single-particle strength function - one-pole approximation 
- is concerned). Consequently, the absolute one-nucleon transfer differential cross section can 
depend, in an important way, on the knowledge of the optical potential inside the nuclear volume. 

The above discussion can also be related to the fluctuation-dissipation theorem (cf. eg. refs. 
iTTSl [741 ). closely related to the Kramers-Kronig dispersion relation. While it is true that single- 
particle motion emerges from the same features leading to collective nuclear vibrations [[6T1l - the 
independent particle shell structure being a result of a collective and concerted motion of all the 
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nucleons allowing free motion to each single of them, and letting themselves felt through the push- 
ing and pulling only when a particle tries to leave the system, forcing it to bounce elastically off the 
nuclear surface - the picture becomes less stable as soon as the particle- vibration coupling strength 
is switched on, leading to real, dressed particles (mechanism also contributing to the imaginary 
part of the optical potential). Within this context one can posit that the associated single-particle 
strength function in general and thus its centroid and width in particular, depend on a delicate 
interplay of spin- and non spin- flip matriK elements, large and small energy denominators, and 
the like. On the other hand the strength function associated with a Cooper pair, being a coherent 
object behaving almost semiclassically, resents much less of the above mentioned inhomogeneous 
damping phenomena. 

Summing up, even without totally consistent optical potentials in the Kramers-Kronig sense, it 
is likely that one can calculate the absolute value of the two-particle transfer cross section between 
members of pairing rotational (vibrational) bands with high accuracy. On the other hand, in the 
case of one-nucleon transfer processes, this possibility is likely to be restricted to the centroid and 
width of the (main peak) single-particle strength function, the integrated area being affected, as 
a rule, by little controllable, non-specific background effects, difficult to estimate and/or remove. 
Such limitations will Ukely be more important in the case of strongly fragmented single-particle 
states. 
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FIG. D.l: Josephson junction between left (L) and right (R) superconductors weakly coupled through a 

thin (DC biased) oxide layer; (b) time dependent nuclear Josephson junction estabUshed at about the 
distance of closest approach by the superfluid nuclei a and A in the reaction a + A — » (a ± n) + A(- + n), 
where n -0, 1,2... is the number of Cooper pairs transferred in the process; (c) dynamical time dependent 
Josephson junction between a closed shell system a displaying strong pairing vibrations (pair addition and 
pair removal mode) and a superfluid nucleus A. In the inset (left) the ground zero point fluctuations of the 
system a associated with the pair addition (arrowed double Une pointing up) and pair removal mode 

(arrowed double line pointing down) are shown. 
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